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ABSTRACT 


This report describes a series of low speed airfoil designs based 
on modifications to the NACA 64-206 airfoil. Designs are based on 
potential flow theory. The report describes one of a series of airfoil 
modifications carried but under Contract NAS 2-8599, Application of 
Multivariable Search Techniques to Optimal Wing Design in Non-Linear 
Flow Fields. Mr. Raymond Hicks of National Aeronautics and Space 
Administration's Aeronautical Division, Ames Research Center, served 
as contract monitor for the present study. 
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AN INVESTIGATION ON THE EFFECT OF 
SECOND-ORDER ADDITIONAL THICKNESS DISTRIBUTIONS TO 
THE UPPER SURFACE OF AN NACA 64-206 AIRFOIL 

by Antony W. Merz and Donald S. Hague 
Aerophysics Research Corporation 

SUMMARY 

An investigation has been conducted on the Lawrence Radiation Center, 
Berkeley, CDC 7600 digital computer to determine the effects of additional 
thickness distributions to the upper surface of an NACA 64-206 airfoil. 
Additional thickness distributions employed were in the form of two second- 
order polynomial arcs which have a specified thickness, y, at a given 
chordwise location, x. The forward arc disappears at the airfoil leading 
edge, the aft arc disappears at the airfoil trailing edge. At the 
juncture of the two arcs, x = x, continuity of slope is maintained. The 
effect of varying the maximum additional thickness and its chordwise 
location on airfoil lift coefficient, pitching moment, and pressure dis- 
tribution was investigated. Results were obtained at a Mach number of 
0.2 with an angle-of-attack of 6° on the basic NACA 64-206 airfoil. All 
calculations employ the full potential flow equations for two dimensional 
flow. The relaxation method of Jameson is employed for solution of the 
potential flow equations. 

Introducing this type of upper surface modification to the NACA 64-206 
airfoil produced results which generally follow trends found previously 
in a similar investigation employing the NACA 64^-212 airfoil. Increases 
in the rearward location of the maximum additional thickness and increases 
in the magnitude of the additional thickness both produce increases in 
the airfoil lift coefficient. Conversely moving the location of maximum 
thickness forward or decreasing the maximum thickness both reduce the 
magnitude of the quarter chord pitching moment. The magnitude of the 
largest pressure peak varies in a complicated manner with maximum additional 
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thickness and its chordwise location. For maximum thickness locations 
forward of the 2/3 chord additional thickness initially produces a re- ; 
duction in pressure peak with a lift coefficient increase. With larger 
amounts of additional thickness pressure peak value and lift coefficient 
tend to rise together. This reversal in C^-Cp trend results from 


max 

the creation of a second peak pressure region aft of the basic airfoil 
leading edge pressure peak. As thickness increases the magnitude of this 
aft peak rises while the magnitude of the leading edge pressure peak 
decreases. Minimum Cp occurs when the two peak pressure values are 


max 

equal. Further increases in thickness beyond this point produce simul- 
taneous increases in both lift coefficient and peak pressure. Solutions 
were difficult to obtain for the thicker airfoils during the present 
study. Increased lift, and its accompanying circulation around the 
relatively sharp leading edge of the NACA airfoil, caused numerical 
difficulties in solution of the potential flow equation for these equations. 
This effect was most pronounced for aft location of maximum additional 
thickness. 

For maximum thickness locations aft of the 2/3 chord location 
additional thickness produces a monotonic rise in both lift coefficient 
and pressure peak magnitude. In these cases the leading edge pressure 
peak always dominates. A consequence of the above behavior is that for 
a given lift coefficient value the peak pressure can be minimized by 
careful selection of the location of maximum thickness and its magnitude. 
Generally as the lift coefficient rises the maximum thickness location 
moves aft. For lift coefficients between 1.0 and 1.6, the chordwise 
location of the maximum thickness varied from 10% to 30%. For this 
increase in lift coefficient, the magnitude of the peak pressure coefficient 
decreases from about 3 to 2.4. 

It should be noted that viscous effects are neglected in the present 
analysis. At the higher lift coefficients the effect of viscosity could 
be significant. Further investigations incorporating a viscous flow model 
are therefore 'desirable. 
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INTRODUCTION 


The National Aeronautics and Space Administration and others are 
currently conducting a series of theoretical and experimental studies 
to define airfoil sections having improved performance from the aspects 
of lift, drag, pitching moment or pressure distribution characteristics, 
references 1 and 2. Analytic investigations using airfoil surface repre- 
sentations based on high-order polynomials may result in impractical 
profiles, for example, very thin trailing edge thickness distributions 
or severe reflexes in the profile. The present study employs low-order 
polynomial arcs of second-order whose characteristics are selected to avoid 
such problems. Previous optimization studies using multivariable search 
techniques, references 1, 3 and 4, generally indicate that shape changes 
which provide increased lift produce unfavorable changes in moment 
characteristics. Conversely profile changes which improve the moment 
characteristics decrease the lift coefficient. With the low-order model 
of the present investigation a systematic examination on the effect of 
profile changes can be carried out. This was accomplished and the trends 
previously revealed by optimization studies were confirmed. An interesting 
by product of the systematic investigation of profile changes and that 
of the previous investigation of the NACA 64^-212 airfoil is that a gain 
in lift coefficient can be produced while reducing the peak negative 
pressures. This tends to decrease the pressure gradient and hence holds 
promise for the development of practical single component high lift 
coefficient airfoils. 

MATHEMATICAL MODELS 

Potential Flow Equation 

Potential flow analysis is based on solution of the two-dimensional 
potential flow equation 

(a 2 -u 2 ) A + (a 2 -v 2 ) 0 - 2uv 0 =0 

( J XX v ' yy xy 

where 0 is the velocity potential, u and v are the velocity components 
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u = 0 X , V = 0 y 

and a is the local speed of sound determined from the energy equation and 
the stagnation speed of sound 

a 2 = a o 2 - (u 2 + v 2 ) 

Solutions are obtained by Jameson’s finite difference scheme, reference 5. 


AIRFOIL PROFILE REPRESENTATION 


Basic Airfoil 


Ordinates for the basic NACA 64-206 airfoil were approximated by 
four cubic chain polynomials in the manner of Hicks 

2 3 

y. = a Fj+ a, x + a 2 x + a^ x ; J - 1,2, 3,4 
J °j A j j j 

Coefficients in the four polynomial arcs are selected on the following 


i = 1 - Arc represents forward portion of upper surface 

F i * V* 

i = 2 - Arc represents aft portion of upper surface 

"j “ 1 . 

i = 3 - Arc represents forward portion of lower surface 


’ ' P 3 * V : 


i = 4 - Arc represents aft portion of lower surface 

H 4 “ 1 

The coefficients a. are determined by introducing four boundary conditions 

on the airfoil profile in each of the four airfoil arcs. Crout's method 
for triangularization and back substitution is used to solve the resulting 
system of linear equations. Note that if four points are specified on the 
aft portion (i* * 2 or 4) , a discontinuity in slope occurs where the poly- 
nomials join. This produces a small ripple in the pressure distribution 
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at the juncture point. However, since the juncture occurs at a region 
of small slope (x = .5) the effect is not significant. The approximate 
NACA 64-206 airfoil developed by this method is presented in Figure 1. 

Additional Thickness 

In the present study additional thickness is limited to the upper 
airfoil surface. The additional thickness has the form 


■ y f. 1 -((ff ] • 

Ay(x) = y £l 1 ; 


x <x 

X> X 


These functions are of second-order varying parabolically with C = |x - x|. 
Additional thickness is zero at the leading edge ( x =0) and trailing edge 
(x = 1) and has a maximum of Ay = y at x = x. Additional thickness and 
slope of the additional thickness are continuous throughout the interval 
0< x<l. The second derivative of the additional thickness distribution 
is constant in the forward and aft airfoil arcs but has a discontinuity 
at the arc junction, x = x. It follows that a continuous polynomial 
representation of the additional thickness distributions, valid in the 
interval 0<x<l would be in the form of an infinite series. This type 
of additional thickness distribution is referred to as a "biquadratic" 
function in recognition of the above characteristics. A sequence of 
biquadratic arcs having varying maximum thickness positions is presented 
in Figure 2. 


PREVIOUS OPTIMIZATION STUDIES 




contour is completely described in terms of the two parameters x and y 
0 = Max [cj = Max £c L (x,y) j 

where 

This two variable multivariable search problem can be solved by a 
combination of directed random-ray and pattern searches, reference 3. 

Table I presents the results of 30 iterations previously obtained in the 
reference 4 study of the NACA 64^212 airfoil using these search procedures 
Lift gains were produced at 27 of the 30 iterations and continued to be 
made at the computation termination. 

Optimization moved the position of maximum thickness to the most 
rearward position allowed, x = 0.9. At termination, lift was increasing 
monotonically with increasing thickness, y. Based on this isolated 
result lift may be maximized for additional thickness of the form assumed 
by moving the position of maximum additional thickness as far aft as 
allowed and introducing as much additional thickness as allowed. 

Moment Minimization 

Minimization of the moment coefficient has the form 

4 . r 

0 = Min [C M ] 

where 

C M = j(x - 1/4) Ap (x) dx 

In previous studies moment minimization resulted in a solution directly 
opposed to lift maximization. The position of maximum thickness moved 
forward and the amount of additional thickness was minimized. Thus the 
basic airfoil tends to have less adverse moment than any airfoil generated 
by addition of biquadratic thickness to the upper surface of the airfoil. 

Optimization Summary 

Optimal airfoil results previously obtained in the reference 4 study 
are summarized in Table II. It can be seen that in all cases previously 
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studied the position of maximum thickness, x, is either at the extreme 
forward or rearward position allowed. Similarly, depending on problem 
specification, the amount of additional thickness should be either 
minimized or maximized. The low dimensionality of this problem (two 
parameters, x and y) permit a ready mapping of results obtained on the 
present modifications to a NACA 64-206 airfoil as a function of x and 
y. This is done in the following section. 

SYSTEMATIC VARIATION OF NACA 64-206 
AIRFOIL SHAPING PARAMETERS 

The present study of a NACA 64-206 airfoil was based on a systematic 
investigation on the effect of variations in the airfoil shaping parameters 
x and y. The resulting airfoils and calculated pressure distributions 
are presented in Figures 3(a) to 3(V). The pressure signatures vary in a 
radical manner with x and y. The basic airfoil exhibits a sharp pressure 
peak at the leading edge. Magnitude of the peak pressure is reduced by 
introducing additional thickness in a forward location, x = .1, and the 
leading edge pressure peak position moves aft. However, if the amount 
of additional thickness is further increased the pressure peak magnitude 
again increases. This peak is located well aft of the leading edge 
pressure peak. This result is due to the creation of a second peak in 
the upper surface pressure distribution. This effect persists until 
rearward locations of x are encountered. For example, (Figure 3(u)), 
introducing additional thickness at x = ,7 theoretically results in a 
rearward "hump” in the pressure distribution somewhat similar to that 
produced by a trailing edge flap. At this extreme aft location the 
increased circulation produced by this pressure hump also produces an 
increased leading edge peak in the airfoil pressure distribution. Flow 
separation would probably be encountered with these rearward additional 
thickness distributions unless devices such as rotating cylinders Or 
blowing were employed. 

Figure 4 illustrates the effect of varying position of maximum 
thickness and maximum thickness on lift coefficient. It can be seen that 
lift coefficient is maximized by increasing both X and y. Generally 
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Figure 4 confirms trends of optimization studies using the NACA 64 1 -212 
airfoil discussed in the previous section. The principal effect of 
changing reference airfoils from the NACA 64^-212 to the NACA 64-206 
is a reduction in the magnitude of pressure coefficients. The dashed 
lines of Figure 4 show the lift coefficient variation of the thicker 
NACA 64^-212 airfoil (reference 4) superimposed on the results obtained 
with the NACA 64-206 airfoil. Lift coefficient values are displaced 
downward by a nearly constant amount when the NACA 64-206 airfoil is 
employed. Since the additional thickness and the basic 6% airfoil thickness 
are additive. Figure 4 presents lift coefficient as a function of thickness. 

To first-order the airfoil thickness required is 

t/c = 6% + y 

As x moves to the extremes of the range the actual airfoil thickness 
is less than this amount as the positions of maximum thickness on the 
basic and additional thickness distributions are significantly different. 

Moment coefficient variation with x and y is presented in Figure 5. 

It can be seen that the increased lift available from additional thickness 
is accompanied by a matching increase in undesirable pitching moment 
coefficient. The conclusion of the previous section that moment coefficient 
is minimized by moving x forward and diminishing y is borne out by Figure 5, 
again confirming previous optimization study results in reference 4. 

As has been noted above, the primary effect of reducing basic airfoil 
thickness from 12% to 6% is to reduce both the lift and the pitching 
moment coefficients. Figure 5 illustrates this effect for pitching moment 
coefficient. The solid line presents results obtained on the NACA 64-206 
airfoil. The dashed lines superimpose previous results obtained with the 
NACA 64 1 ~212 airfoil. 

Figure 6 provides another means of studying the simultaneous variations 

of C. and C . That is, the desirable characteristics of high lift and 
L M 

low moment are attainable only in a relative or weighted sense (reference 4) . 

For the NACA 64-206 airfoil, the C T -C w variation at a constant value of 

L M 

x is only nearly linear. The desired slope of this line is as small as 
possible. As shown in Figure 6, this is attained by making the point 
of maximum additional thickness close to the leading edge. The numeral 
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procedure, however, may not converge when the additional thickness is 
added too close to the leading edge. The dashed line of Figure 6 
corresponds to x = .2, and the data extend below this line for x = .1 
only at y = .03 and .06 due to these numerical difficulties. However, 
the trend is clear when the C^-C^ trade-off is measured by the line function 
criteria of reference 4. Maximum thickness should be employed and intro- 
duced as near to the leading edge as possible. 

Figure 7 presents the relationship between pressure peak and lift 

coefficient fc * a range of x and y values. For each value of y (maximum 

additional thickness) there is a point at which the pressure peak magnitude 

is minimized. Cross plotting the peak pressures as a function of C in 

L 

Figure 8 reveals the minimum peak pressures as a function of C . 

L 

Figure 9 plots the position of maximum additional thickness as a 
function of C^. As increases, x moves aft. The associated values of 
y required for the low peak pressure is also plotted in Figure 9. The 
amount of additional thickness required for a minimum pressure peak 
increases with C^. Finally, Figure 10 plots the minimum Cp attainable 
as a function of C L using the biquadratic additional thickness airfoil 
model. For the range of thicknesses studied here, Cp, reduces with 

max 

C L flattening out at C^ = 1.5. Minimum C p values are higher than those 
attained in the Reference 4 study using the 64^-212 airfoil. 

CONCLUSION 

A numerical investigation into a class of modified airfoil shapes 
has been completed using full two-dimensional flow potential flow equations. 
Airfoils studied were obtained by modifying the NACA 64-206 airfoil by 
additional thickness distributions based on a biquadratic variation with 
chordwise position. Free stream Mach number was held constant at M = 0.2 
and the basic airfoil is held at 6° angle-of-attack. Results of the 
study may be summarized as follows and are generally in agreement with 
reference 4: 

1. Significant changes in pressure distribution, lift and 
pitching moment can be introduced by the biquadratic 
thickness modification. 
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The requirements for improving lift and moment coefficient 
characteristics are directly opposed tb each other. That 
is, increases in lift result in increases in adverse moment. 
Conversely, decreases in adverse moment produce decreases in 
lift. 

High lift airfoils require the addition of a thickness 
distribution biased to the rear of the foil and as much 
thickness addition as possible. 

Low adverse moments require a thickness distribution biased to 
the front of the foil and as little additional thickness as 
possible. Therefore, the best airfoil based on moment con- 
siderations is the unmodified foil. 

Favorable lift/moment trade-off characteristics are obtained 
by a thickness distribution biased to the front of the foil 
employing as much thickness as possible. 

There exists a class of airfoil exhibiting low peak pressures 
for a given C which require an intermediate location of 

Li 

maximum additional thickness and thickness amount. Generally, 
the position of maximum additional thickness moves to the 
rear with increasing C^, and the amount of additional thickness 
required increases with increasing C^. •" 
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FIGURE 3(k) . MODIFIED 64-206 AIRFOIL, x = .4, y = .12 



FIGURE 3(1). MODIFIED 64-206 AIRFOIL, x = .5, y = .06 



FIGURE 3(m). MODIFIED 64-206 AIRFOIL, x = .5, y = .09 
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TABLE I 


CONVERGENCE FOR C T MAXIMIZATION 
L 


M 

C 

JJJ 

0 

1 

1 

1° 

1 

2 

10 

1 

3 

10 

1 

5 

2 

I 

6 

2 

1 

7 

2 

i 

3 

2 

1 

9 

1° 

2 

10 

10 
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11 

10 

2 

12 

2 
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13 

2 
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2 

2 

15 

w 2 

2 

16- 

10 

3 

17 

IP 

3 

16 

10 

3 

20 

2 

3 

21 

2 

-3- 

22 

2 

3 

23 

2 

3 

24 

!o 

4 

25 

To 

4 

26 

To 

n 

“l 

2 7 

U 1 

2 

4 

26 

2 

4 

29 

2 

4 

30 


X 


ALPHA ( 1) 
.2500 
. 2506 
.2583 
,2626 
,2753 
,2879 
.3132 
,3638 
#36 88 
. 3758 
.3824 
.4010 
.4196 
.4567"“ 
.5311' 
.5505 
,5675 
,5808 
. 6305 
* 680 2 . . 
,7796 
,9000 
,9006 
,9000 
.9000 
.9000 
.9000 
,9000 


y 


ALPHAf 2) 
\‘„ 0000E-03 
1 , 0365E-03 
T. J 667E-03 
r.4 155E-0 3 
l’.8311E-03 
2, 2466E-03 
3'.O776E-03 
4.7398E-03 
4'. 7 53 3E- 0 3 
5‘. 2172E-03 
5,561 IE-03 
6'. 3824E-03 
7 . 2 0 3 7 E - 0 3 
8'.3462E-03 
1 .2131E-02 
1 .2761E-02 
1 .2B12E-02 
1.3494E-02 
l', 4656E-02 
T.6219E-02 
1 , 09 4 4 E- 02 
2 . 4 394 E- 0 2 
2 5 0 1 6 E " 0 2 
2, 7157E-02 
3.01 16E-02 
3.5838E-02 
4 , 1 56 0 E "02 
5*. 3003E-02 


FUNCTNf 1 ) 
-‘.8941 
-.8942 
-.5952 
-.8968 
-,8982 
-.9022 
-.9078 
-.9196 ' 
-.9196 
-.9216 
-.9258 
-.9320 
-.9385 

-.9520 

-.9826 

-.9863 “ 

-.9920 

-.9993 

-1.018 

-1.041 

-1.105 

-1.303 

-1.313 

-1,350 

-1.399 

-1.493 “ 

-T.586 

-1,765 


